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The existing experimental data on total cross sections for hadron production in e + e~ annihilation 
in the resonance region yfs = 3.8 - 4.8 GeV, usually presented in terms of the parameter R, are 
critically examined. It is shown that the Crystal Ball and BES measurements are in excellent agree- 
ment, and their analysis leads to consistent resonance parameters for the three vector resonances 
above the DD threshold. The results for the widths are found to be considerably different from the 
presently adopted values, and have much smaller errors. 
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PACS numbers: 13.65.+i, 13.20.Gd, 13.25.Gv 



One of the most important observables in particle 
physics is the total cross section for hadron production 
in e + e~ annihilation, a(e + e~ — > hadrons). It is gener- 
ally presented as the ratio to the QED muon production 
cross section a(e + e~ — > muons). The ratio is conceptu- 
ally very simple for the continuum cross sections. In the 
quark model, 
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with the sum extending over the quark flavors open 
upto the cm. energy y/s. In QCD the ratio acquires 
corrections in powers of (a s /7r). In the limit of massless 
quarks, the corrections are known to the third order, and 
the expansion is satisfactorily convergent. 
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For four open flavors (u, d, s, c), C% = 1.5243, and C3 = 
-11.520 [1,2]. Because of this well established correction 
factor, measurements of R(\/s) are often used to deter- 
mine the strong coupling constant a s . Above the open 
flavor thresholds, measurements of R are used to deter- 
mine parameters of vector resonances. In particular, the 
totality of our present knowledge about charmonium res- 
onances above the DD threshold at 3739 MeV comes 
from a single measurement of R(y/s) by DASP [3]. This 
paper is devoted to an alternate analysis of the presently 
available R(\/s) data for the parameters of these reso- 
nances. 

The experimental determination of R generally con- 
sists of measuring the total cross section for the produc- 
tion of hadrons, and dividing it by the lowest order QED 
cross section, a(e + e~ — > muons) (in nb) = 47rag m /3s = 
86.8/s (GeV 2 ). 

The existing data for R, which span the region y/s = 
3.8 - 4.8 GeV, which contains the resonance structures 



TABLE I: R Measurements from different sources [Ref. 3-7]. 



Ref. 

(Year) 


(GeV) 


Steps 
(MeV) 


%Stat. 
Error 


%Syst. 
Error 


DASP [3] 
(1978,9) 


3.6-4.8 


10-15 


5-10 


~ 12 


MARK I [4] 


3.4-4.7 


10 


5-30 


20-10 


(1982) 










PLUTO [5] 


3.6-5.0 


20-40 


3-9 


10-15 


(1982) 










GB[6] 


3.6-4.5 


6 


3-7 


~ 9 


(1986) 










BES [7] 


3.7-4.8 


5-10 


~ 4 


~ 5 


(2002) 











TABLE II: Parameters of the higher -0(1 ) resonances. 



Source 


Mass 


Ttot 


Tee 




(MeV) 


(MeV) 


(keV) 


DASP [3] 


4040 ± 10 


52 ± 10 


0.75 ±0.15 


(1978 


4159 ± 20 


78 ±20 


0.77 ±0.23 


-1979) 


4417 ±20 


66 ± 15 


0.49 ±0.13 


MARK I [9] 


4417 ± 10 


33 ± 10 


0.44 ±0.14 


(1976) 








PDG[1,10] 


4415 ± 6 


43 ± 15 


0.47 ±0.10 


(1980-2002) 
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of interest, are listed in Table I along with their princi- 
pal characteristics [3-7]. Radiative corrections have been 
made in all data sets and heavy lepton contributions have 
been removed. The numerical data for DASP [3] measure- 
ments have been obtained from Ref. [8]. 

Because we wish to determine Breit-Wigner resonance 
parameters we have converted the R values to cj(e + e _ ) in 
nb = 86.8 R/s (in GeV). These cross sections are plotted 
in the four panels of Fig. 1. We do not use the PLUTO [5] 
data because the rather coarse spacing of its data points, 
20-40 MeV, does not allow analysis for resonances whose 
widths arc of the same order. Only DASP [3] has re- 
ported a resonance analysis of their data, and in Fig. 1 
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FIG. 1: Total hadron production cross sections in nb obtained 



as a(e^ 



hadrons) = 86.8 R/s (GeV) from the R mea- 



surements of DASP [3], MARK I [4], Crystal Ball (CB)[6], and 
BES[7]. Only statistical errors are shown. The common curve 
superposed on all four data sets is that obtained by DASP [3] 
as the best fit to their data. 

we have superposed on all data sets the curve which rep- 
resents the fit obtained by DASP [3] by fitting their data 
with three Breit-Wigner resonances and a background 
smoothly connecting their measured R values at the two 
extremes of 3.6 GeV and 5.0 GeV. The \ 2 /d.o.f. val- 
ues listed in the panels of Fig. 1 indicate that the fits 
with DASP parameters are quite poor in all other cases, 
particularly to the Mark I data. The resonance param- 
eters obtained by DASP are listed in Table IF The Par- 
ticle Data Group (PDG) adopted the DASP parameters 
for the resonances at 4.04 and 4.16 GeV, and averaged 
the DASP parameters for the third resonance with those 
given in an early report by MARK I [9]. These are the 
only available resonance parameters, and they have re- 
mained unchanged from PDG-1980[10] to PDG-2004[1]. 

The three pre-1990 data sets shown in Fig. 1 differ 
both qualitatively and quantitatively. This was noted 
early by, for example, Eichtcn et al. [11], and seems 
to have discouraged theoretical attempts to understand 
higher vector states of charmonium. One of the most ob- 
vious differences is that neither MARK I, nor CB, nor 
BES find the deep minimum at ~ 4.1 GeV observed in 
the DASP data[12]. Similarly, all three appear to observe 
broader structures than those claimed by DASP, and in 
the initial report by MARK I [9]. There are also obvious 
differences in the absolute levels of cross-sections which 
can be best compared in the off-resonance regions. In 
the region = 3.60 — 3.65, Mark I [4] cross sections are 
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FIG. 2: Total hadron production cross sections and the 
X 2 /d.o.f for the best fits shown: (Top) Crystal Ball (CB)[6], 
(Bottom) BES [7]. 

~ 27% larger than those of DASP [3]. Similarly, in the 
off- resonance region, 6.0—7.5 GeV, where additional data 
from CB are available [13], Mark I data [4] are <~ 28% 
larger than CB. The only two data sets which arc not 
only in qualitative, but excellent quantitative agreement, 
are those due to CB [6] and BES [7]. They also have 
the smallest step size and the smallest statistical and sys- 
tematic errors. We base the present resonance parameter 
determination on these two data sets. 

In order to obtain the best resonance parameters pos- 
sible, we fit the separate CB and BES data sets, each 
with three Breit-Wigner resonances, with 

o-bw{V~s) = (3tt /Ap 2 )T ee r h /[(^7s - Mf + (r 4ot /2) 2 ], 

where T ee , Th, and T to t are the mass independent elec- 
tron, hadron and total widths, respectively, for a vector 
resonance of mass M produced in the head-on collision 
of e + and e~ , each of momentum p. There is little ev- 
idence in the data for any substantial variation of the 
continuum background in this region of i/s. We there- 
fore parametrize it with a linear function, a(bkgd) = 
A — B(y/s — 3.8 GeV). The results of the separate best fits 
are listed in Table III and the fits are shown in Fig. 2. 
As expected, the possible normalization differences be- 
tween CB and BES, and their individual absolute cross 
section uncertainties, quoted as systematic errors in Ta- 
ble 1, do not reflect in the masses or the total widths; 
they are found to be in near perfect agreement. Also, 
the differences in electron widths are well within errors. 
It is very gratifying that two independent measurements, 
made twenty years apart, are in such good agreement. 
Because of this excellent agreement in the results of the 



3 



TABLE III: Summary of results. Masses M (i) and total 
widths T[ % J t are in MeV, electron widths T^} are in keV, the 
parameter A is in nb, and the parameter B is in nb/GeV. The 
final row shows the xVd.o.f. for each fit. 





M w 


r (l) 

1 tot 


r (i) 




(MeV) 


(MeV) 


(keV) 


PDG[1] 


4040 ± 10 


52 + 10 


0.75 + 0.15 


CB[6] 


4037 ± 2 


85 + 10 


0.88 + 0.11 


BES[7] 


4040 ± 1 


89 + 6 


0.91 + 0.13 


CB+BES 


4039.4 ± 0.9 


88 + 5 


0.89 ± 0.08 




M< 2 > 


r (2) 
L tot 


r (2) 

J- ee 


PDG[1] 


4159 + 20 


78 + 20 


0.77 + 0.23 


CB[6] 


4151 + 4 


107+10 


0.83 + 0.08 


BES[7] 


4155 ± 5 


107+16 


0.84 + 0.13 


CB+BES 


4153 ± 3 


107 + 8 


0.83 + 0.07 






r (3) 
1 tot 


r (3) 
-L ee 


PDG[1] 


4415 ± 6 


43 + 15 


0.47 + 0.10 


CB[6] 


4425 ± 6 


119 + 16 


0.72 + 0.11 


BES[7] 


4429 ± 9 


118 + 35 


0.64 + 0.23 


CB+BES 


4426 ± 5 


119 + 15 


0.71 + 0.10 




A 


B 


x 2 


DASP[3] 
CB[6] 


Polynomial 
14.2 + 3.5 


1.5 + 0.4 


2.1 
0.99 


BES[7] 


13.7 + 4.5 


1.5 + 0.5 


1.23 



CB and BES measurements, we consider their weighted 
averages, given in the last column of Table III, as our 
final results. 

From Table III it is apparent that the presently deter- 
mined masses of the three resonances have much smaller 
errors than the masses determined by DASP [3] , but are in 
general agreement with them. However, the total widths 
determined in the present analyses are quite different 
from those determined by DASP, and those adopted by 
PDG2004. Our values of r«, and and - 67%, 
37%, and 179% larger, respectively than the correspond- 
ing PDG values. The corresponding electron widths de- 
termined by us are also larger, by 23%, 8%, and 51%. 
The errors in our width determinations are also more 
than factor two smaller. It is also worth noting that the 
present analysis does not rule out the presence of narrow 
structures other than the three broad ones analyzed here. 

It is interesting to note that the average value of R over 
the entire region of measurements is < R >= 3.63 ± 0.13 
for CB, and < R >= 3.35 ± 0.14 for BES. These values 
are close to the prediction, R = 3.67 obtained from eq. 



2 for four flavors, with a s = 0.30, which corresponds to 
a s (m z ) = 0.117[1]. 

Eichten et al. have identified, and most later poten- 
tial model calculations agree with their identification, 
that the three resonances are 3 3 Si (4039 MeV), 2 3 L>i 
(4153 MeV), and 4 3 S*i (4426 MeV). Three early poten- 
tial model calculations have made predictions of the lep- 
tonic widths of 3 3 Si and 4 3 Si states. Quigg and Rosncr 
predicted r ee (4039) = 1.00 keV, and T ee (4426) = 0.51 
keV for their logarithmic potential [14]. Eichten et al.[ll] 
predict r ee (4039) = 1.5 keV, and r ee (4426) = 1.1 keV 
with the Cornell potential, and Buchmuller and Tye[15] 
predict r ee (4039) = 1.68 keV, and r ee (4426) = 1.31 keV 
(assuming T ee (J/xjj) = 5.25 keV) with a QCD based po- 
tential. In principle, all other potential model calcula- 
tions give predictions for wave functions and therefore 
leptonic widths, but no subsequent explicit predictions 
can be found in the published literature. 

Theoretical studies of the hadronic properties of the re- 
gion above the DD threshold are extremely rare, mainly 
due to the fact that there were "obvious disagreements 
between experiments" [11]. Despite this problem, Eichten 
et al. [11] made a rather detailed, though unavoidably 
conjectural, study of this region. They attempted to ex- 
plain the nature of the observed structures in the 3.9 to 
4.2 GeV region in terms of the 3S and 2D resonances, the 
mixing between them, and the opening of the DD* and 
D*D* thresholds at 3880 MeV and 4020 MeV respec- 
tively. They did not succeed in reproducing the deep 
minimum observed in the data of DASP, but as we have 
seen, the CB and BES data definitely rule out the min- 
imum, and our resonance parameters fit these data very 
well. Eichten et al. explained the structure at 4.4 GeV 
entirely in terms of the AS resonance and its D*D* de- 
cay. The discovery of strange D s and D* came later, 
and the opening of the D S D S and D* S D* S thresholds at 
3940 MeV and 4224 MeV, respectively, can be certainly 
expected to play a significant role in the understanding 
of the observed structures, particularly that of the 4S 
state at 4426 MeV. We hope that our clarification of the 
experiment data, and the more precise determination of 
the resonance parameters will lead to renewed theoretical 
interest in the understanding of the challenging physics 
of these higher vector states of charmonium. 
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